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We report the discovery of WASP-13b, a low-mass M p = 0.46 '_nm^J transiting exoplanet with an orbital period of 4.35298 ± 0.00004 days. The 
transit has a depth of 9 mmag, and although our follow-up photometry does not allow us to constrain the impact parameter well (0 < b < 0.46), 
with radius in the range R p ~ 1.06 - 1.21/fy the location of WASP- 13b in the mass-radius plane is nevertheless consistent with H/He-dominated, 
irradiated, low core mass and core-free theoretical models. The G1V host star is similar to the Sun in mass (M, = 1.03*q^M o ) and metallicity 



([M/H]=0.0 ± 0.2), but is possibly older (8.5^ Gyr). 
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The discovery of transiting planets is a prominent theme in mod- 
ern astroph ysics. Even in the era of space-borne surveys such 
as CoRoT dBarge et al. 2 007). the detection of a new transiting 
planet remains an important and celebrated discovery. In part 
this is because these are the only systems for which accurate 
physical parameters can be determined, which in turn enables 
their mass-radius relationship to be used as a diagnostic to con- 
strain models of exoplanet structure and evolution. Currently 
there are some sixty transiting exoplanets known, and these man- 
ifest significant diversity in their physical properties. 

There are four leading ground-based transit surveys: 
HATNet (jBakos et al. 2004 . the Trans- Atlantic Exoplanet 
Survey (iDunham et al. 2004 lO'Donovan et al. 2006b. WASP 
dPollacco et al. 20061) and XO (iMcCullough et al. 2006b . Each 
uses specialist instruments capable of imaging the sky over ex- 
tremely large angular scales, and consequently are optimised 
to obtain high precision photometry for relatively bright stars. 
Early p redictions w ere optimistic about the expected planetary 
yields (Home 2003) in such wide-angle surveys. As a clearer 
understanding of th e effects of systematic noise on th e photom- 
etry has emerged dPont et al. 20061: ISmith et al. 2006b . the dis- 
covery rates now broadly reflect these predictions. For exam- 
ple, the WASP survey announced the discovery of thirteen tran- 
siting exoplanets over the period August 2007 to April 2008. 
While space-based surveys are superior to ground-based ones for 
the detection of small planets and long-period systems, ground- 
based surveys are likely to remain important as their extremely 
large fields-of-view make them ideal for surveys of bright stars, 
which are well suited for detailed follow-up observations using 
other facilities. 

In this paper we describe the discovery of a new, relatively 
low-mass exoplanet, WASP- 13b, which was detected as part of 
the SuperWASP-North survey. 



2. Observations and Data Reduction 

The WASP Cameras are wide-field imaging facilities designed 
for the detection of exoplanetary transits. There are two simi- 
lar facilities: SuperWASP-North (hereafter SuperWASP-N) on 
the island of La Palma in the Canary Islands, and WASP- 
South located at Sutherland, South Africa. The instrumen- 
tation and infrastructure used to obtain, store and reduce 
the data are described in detail in Pollacco et al. (2006). 
WASP-13 (= IS WASP J092024.70+335256.6 = 2MASS 
J0920247 1+3352567 = USNO-B1.0 1238-0183620), a V = 
10.51 G1V star in Lynx, was monitored with the SuperWASP- 
N Camera from 2006 November 27 to 2007 April 1, during 
which 3329 30-second images were obtained with a cadence 
of ~ 7 minutes. It was identified as a transit candidate using 
the algorithm outlined by Collier Cameron et al. (2007), and its 
lightcurve is shown in Figure[T](top panel). 

Further /?-band photometry was obtained with the 0.95 m 
James Gregory Telescope (JGT) located at St Andrews, 
Scotland, during the transit of 2008 February 16. The camera 
on this telescope consists of a 1024x1024 e2v CCD with an 
unvignetted field-of-view of 15' in diameter. A total of 1047 
15-second exposures was obtained in clear conditions with see- 
ing of 4 - 5" throughout the night. The last 20 images in 
the sequence were obtained under thick cloud cover and are 
not used in the final lightcurve. Data were processed using 
the Cambridge Astr onomical Survey Uni t data reduction and 
photometry pipeline dlrwin & Lewis 200lh . Aperture photome- 
try with an 8" radius aperture was performed on all stars in the 
field, and seven bright reference stars were selected for the dif- 
ferential photometry. The flux of the combined reference star is 
dominated by a single bright star, HD80408, situated 7.6' from 
the target and approximately two magnitudes brighter. HD80408 
(V-K=1.9) is within one spectral class of the planet host star (V- 
K=1.4). 



2 



I. Skillen et al.: The 0.5Mj transiting exoplanet WASP- 13b 



Table 1. Log of SOPHIE observations of WASP-13. 



BJD 


Phase 


SN26 


RV 




v 

v span 








(km/s) 


(km/s) 


(km/s) 


2454508.48441 


0.875 


45.7 


9.8660 


0.0088 


0.020 


2454509.45978 


0.099 


56.1 


9.8126 


0.0072 


-0.005 


2454510.46914 


0.331 


51.7 


9.8022 


0.0078 


0.017 


2454510.55389 


0.350 


45.3 


9.7634 


0.0088 


0.007 


2454511.32838 


0.538 


51.1 


9.8365 


0.0080 


0.016 


2454511.55075 


0.579 


50.1 


9.8690 


0.0082 


0.009 


2454511.57290 


0.584 


47.3 


9.8492 


0.0086 


0.032 


2454511.64730 


0.602 


46.6 


9.8667 


0.0088 


-0.003 


2454512.42051 


0.779 


60.9 


9.9062 


0.0068 


0.028 


2454512.55184 


0.809 


51.2 


9.8770 


0.0080 


0.010 


2454512.65565 


0.833 


45.0 


9.8829 


0.0090 


-0.013 



The transit ingress was adversely affected by extinction be- 
cause the data were obtained at high airmass in humid con- 
ditions. The data display a trend with airmass which we re- 
moved by subtracting a third-degree polynomial fit to the resid- 
uals from an initial model fit to the original differential light 
curve. There was essentially no correction to the measurements 
obtained at secZ < 1.5, and the amplitude of the correction for 
1.5 < secZ < 1.8 ranges from 0.5 to 5 mmag. The adopted JGT 
lightcurve is shown in Figure Q] (middle panel). 

WASP-13 was observed with the Observatoire de Haute- 
Prove nce 1 .93 m telescop e and the SOPHIE stabilised spectro- 
graph dBouchv et al. 2006h during the period 2008 February 11- 
15, and eleven spectra were acquired with exposure time of 600- 
seconds and signal-to-noise ~ 40 - 60 measured in echelle or- 
der 26 at A ~ 550 nm (SN26). We configured the instrument in 
its high efficiency mode with a resolution of R =40000, acquir- 
ing simultaneous star and sky spectra through separate fibres. 
Thorium- Argon calibration images were taken at the start and 
end of each night, and at 2- to 3-hourly intervals throughout each 
night. The radial-velocity drift never exceeded 2-3 m/s, even on 
a night-to-night basis. Spectra were reduced with the standard 
SOPHIE pipeline and corrections applied for lunar contamina- 
tion of the cross-correlation function as needed. 

A further three spectra of WASP-13 with resolution 
R=46000 were acquired with the 2.5 m Nordic Optical 
Telescope (NOT) using the FIES echelle spectrograph, with the 
aim of determining the host-star parameters (see Section 13.11 ). 
These spectra were extracted with the bespoke data reduction 
package, FIEStool. 



3. Results 

The SuperWASP-N and JGT lightcurves show the presence of a 
9 mmag dip of duration ~ 3.9 hours which repeats with a period 
of ~ 4.35 days. 

The radial velocities derived from the SOPHIE observations 
are listed in Table [1] and plotted in the lower panel of Figure Q] 
WASP-13 exhibits radial-velocity variability in phase with that 
expected from reflex motion caused by a transiting exoplanet. 
We examined the line-bisector span, V spa n (Table 1), in the man- 
ner described by Christian et al. (2009) to search for asymme- 
tries in spectral line profiles that could result from unresolved 
binarity or indeed stellar activity. Such effects would cause the 
bisector spans to vary in phase with radial velocity, but no sig- 
nificant correlation is detected. We conclude that the observed 
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Fig. 1. Phased light and radial-velocity curves for WASP-13 ob- 
tained with SuperWASP-N (top), the JGT (middle) and SOPHIE 
(bottom). MCMC model solutions are shown as dotted lines 
for each data set, and the radial velocity model includes the 
predicted Rossiter-McLaughlin effect assuming the upper-limit 
vsin/ (4.9 km/s) 



photometric and radial-velocity variability is caused by an orbit- 
ing, planet-mass body. 

3.1. Stellar and Planetary Parameters 

We used the NOT echelle spectra to derive the host-star param- 
eters (Table |2). Following proced ures developed in ou r analyses 
of similar systems e.g. WASP-1 (Stempels et al. 20071) . we find 
r eff = 5826+ 100 K, log g = 4.04+0.2 dex and [M/H] = 0.0+0.2 
dex, which is consistent with solar metallicity. The effective tem- 
perature and near-solar mass (below) suggest a Gl spectral type. 
The host star has a detectable Li 6708 line from which we derive 
a lithium abundance log n(Li) = 2.06 + 0.1 dex. The rotational 
line profile of WASP-13 is unresolved in the FIES spectra, and 
we derive an upper limit to v sin i of 4.9 km/s by subtracting in 
quadrature the macroturbulence appropriate to T e ff = 5826 K 
(4. 1 km/s) from the instrumental profile (6.4 km/s). 

In addition to the spectrum analysis we used photome- 
try from Tycho, V T = 10.51 and (B - V) T = 0.89, and 
2MASS, (V T - H) = 1.33 and (V T - K) = 1.39, to esti- 
mate the effective tempera ture using the Infrared Flux Method 
dBlackwell & Shallis 19771) . This yields 7/ eff = 5935 + 183 K, in 
close agreement with that obtained from the spectroscopic anal- 
ysis. The Tycho and 2MASS colo urs suggest a spectral type of 
F9V (ICollier Cameron et al. 20061) . 

The light and radial-velocity curves for WASP- 1 3 were mod- 
elled simultaneously using the method described by Pollacco et 
al. (2008). The initial solution from the Monte-Carlo Markov 
Chain (MCMC) routine converged with a stellar density p» = 
0.43^'jy p Q . To determine the mass and age of WASP13 we 
compared its structure and effective temperature with the solar- 
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Table 2. Parameters of WASP-13. 



R.A. = 09*20 m 24 s .70, Dec. = +33°52'56".6 



Teg 

[M/H] 
log n(Li) 
v sin i 

Spectral type 
V mag 
Distance 



5826 ± 100 K 
4.04 + 0.2 
0.0 ± 0.2 
2.06 ±0.1 
< 4.9 km/s 
G1V 
10.42 
155 ± 18 pc 



metallicity stellar evolutionary models of Girardi et al. (2000). 
In Figure 2 we plot the inverse cube root of the stellar density 

— 1/3 1/3 

p* — Rt/MJ (solar units) against effective temperature for 
the model mass tracks and isochrones, and for WASP13. We 

— 1/3 

adopt this parameter space because p* unlike R* or luminosity, 
is measured directly from the light-curve and is independent of 
the effective temperature determined from the spectrum (Hebb et 
al. 2009). We interpolated the evolutionary tracks and isochrones 

— 1/3 

in the p, - T e g plane and find the mass of WASP 13 to be 
M, = 1.03+o ( 'gM and its age to be 8.5^| Gyr. Uncertainties 
in the derived stellar density, temperature and metallicity are in- 
cluded in the overall errors on the age and mass, but systematic 
errors due to differences between various evolutionary models 
are not. The large error in metallicity of +0.2 dex contributes 
significantly to the uncertainty in the mass and age, and a more 
accurate spectral synthesis would improve the precision of these 
parameters. Nevertheless, we re-ran the transit fitting code a sec- 
ond and final time, adopting an initial value for the stellar mass 
of 1.03 M Q and assuming a 10% uncertainty in this parameter. 
Our results are summarised in Table [3] 

The adopted solution yields an impact parameter b = 
0.46^2j, orbital inclination i = 86.9^ 2 degrees and transit 
duration r = 0.163 + 0.003 day, leading to a radius of R» = 
1.34+j}[j4/? for the 1.03 M host star. According to the stellar 
models, a solar metallicity star of this size and mass has evolved 
off the zero-age main sequence and is in the shell hydrogen burn- 
ing phase of evolution with an age of 8.5 Gyr. The Li abundance 
measured in the spectral synthesis also suggests the star is sev- 
eral Gyr old, but it does not provide a precise age determination. 
The abundance is similar to, or slig htly less than, levels fou nd in 
open clusters with ages of 2-8 Gyr (Ses tito & Randich 2005b . 

Our data suggests a large stellar radius and an old age for 
the host star. However, the JGT photometry does not constrain 
the impact parameter strongly. This affects the derived stellar ra- 
dius and as a consequence, the derived age and planetary radius. 
Therefore, we also present a solution for the b — case, which 
gives a lower limit to the stellar and planetary radii. We note that 
the x 2 with respect to the b = model fit is only marginally 
worse than the overall best fitting model. We encourage acquisi- 
tion of higher quality photometry of this object to enable more 
accurate host-star and planet radii to be determined. 

Although we provide only relatively weak constraints on the 
planet's radius, its mass is well constrained from the radial- 



velocity analysis. With M p 



0.46--M y , 



WASP- 13b is 



amongst the lowest-mass transiting exoplanets, and with an in- 
flated radius in the range R p ~ 1.06—1.212?/ its position 
in the mass-radius plane is broadly consistent with the H/He- 
dominated, low core mass and core-free irradiated models of 
Fortney et al. (2007). 
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Fig. 2. Location of WASP-13 in the (p/p Q r I/3 vs T eff (K) plane 
compared to solar-metallicity stellar evolution mass tracks from 
Girardi et al. (2000). The mass tracks are labelled, and the 
isochrones are 0. 1 Gyr, solid; 1 Gyr, dashed; 5 Gyr, dot-dashed; 
10 Gyr, dotted. According to these models WASP-13 has a mass 
of M, = 1.03!° } 'M and an age of 8.5 t^Gyr. 



The Safronov number for WASP- 13b is = 0.039 ± 0.008 
for the b = 0.46 case, and = 0.043 + 0.007 for the b = 
case. This places it amongst the hotter, less-massive Class II gi- 
ant exoplanets in the classification scheme of Hansen & Barman 
(2007), who proposed that transiting exoplanets can be distin- 
guised into two classes according to their equilibrium tempera- 
tures and Safronov numbers. Recently, Torres et al. (2008) re- 
ported additional support of this dichotomy. However, Fressin et 
al. (2009) suggested on the basis of simulations of a model pop- 
ulation of stars and exoplanets that the apparent grouping into 
distinct classes is not statistically significant. A bimodal distri- 
bution of Safronov number for giant exoplanets would have im- 
plications for models of their formation and evolutionary history. 
Continued discovery and parametrization of new transiting exo- 
planets is needed to confirm this hypothesis. 
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